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Abstract 

We present the asymptotically fastest known algorithms for some basic problems 
on univariate polynomial matrices: rank, nullspace, determinant, generic inverse, 
reduced form (HI E3 EES! El- We show that they essentially can be reduced to two 
computer algebra techniques, minimal basis computations and matrix fraction ex- 
pansion/reconstruction, and to polynomial matrix multiplication. Such reductions 
eventually imply that all these problems can be solved in about the same amount of 
time as polynomial matrix multiplication. 
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00 1 1 Introduction 

O ' 

We aim at drawing attention to today's asymptotically fastest known algorithms for com- 
puting with polynomial matrices. In particular, we shall focus on the following problems: 
compute the rank, a right or left nullspace, the determinant, the inverse and a column- or 
row-reduced form of a given polynomial matrix. Polynomial matrices are quite common in 
^ ■ the analysis of multivariable linear systems and Kailath's treatise Linear Systems [TU] is a 
good illustration of this. 

Recently, algorithms have been designed |Hl El UE1 EEI! that allow to compute solutions to 
these problems in essentially the same amount of time as when multiplying two polynomial 
matrices together. More precisely, given a field K — for example the complex numbers, the 
rationals or a finite field — and given a polynomial matrix A e K[x] nxn whose entries have 
degree in x bounded by d, these algorithms allow to compute rank A, ker/L, det A and 
to row-reduce A in CT^n^d) operations in K, and to compute A -1 when A is generic in 
(J~(n 3 d) operations in K. Here, CT^^d) is the best known asymptotic bound for multiplying 
two matrices in K[x] nxn of degree d El, where 2 < uj < 2.376 is the exponent of 
matrix multiplication over K jH Chapter 15]. Using schoolbook matrix multiplication, we 
have lo = 3 and the bound 0~(n^o?) becomes (T(n 3 d). Furthermore, the soft-0 notation 
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0~ simply indicates some missing logarithmic factors of the form a(logn) (logn) 7 for 
three positive real numbers a,/3,j. By achieving the complexity estimate 0~{n u) d) ) these 
algorithms improve upon all the complexity estimates that were known previously 
In this paper, evidence is given that the key tools for such improvements are: 

• Minimal bases of K[x] -modules; 

• Expansion/reconstruction of polynomial matrix fractions. 

The former has the same flavour as in [H] while for the fractions we heavily rely on the 
concepts in ^01 Chapter 6]. Two kinds of minimal bases, namely approximant bases and 
nullspace bases, are studied in Sectional There we will see that such bases are small enough 
to be computed fast, that is, in 0~(n w d) operations in K. Polynomial matrix fractions are 
matrices F G K(x) nxn , where K(x) is the field of rational functions over K. By expansion of 
F, we thus mean a power series expansion F = ^ iX% ^ K[[x]] nx " and by reconstruction 
of F we mean a left or right quotient of polynomial matrices like F = A~ 1 B or F = BA^ 1 . 
It turns out that all we need is truncated expansions and reconstructed quotients of rather 
low degree, both of which can be computed fast as seen in Section 01 The key idea here is 
that an approximant of sufficiently high order — with respect to the input problem — may 
lead to an exact solution over K[x}. This is well-known in computer algebra, at least for 
scalar rational functions §5.7], but as far we know the extension to the matrix case is 
more recent [HI El EEE1 Ej • 

Minimal bases and matrix fractions are interesting not only because they can be com- 
puted fast, but also — and, perhaps, mainly — because computing a minimal basis and ex- 
panding/reconstructing a matrix fraction are problems to which we can reduce all other 
problems like rank, left nullspace, determinant, generic inverse and row-reduced form. The 
goal of Section 0] is precisely to show this: there the above problems are thus seen as 
applications of the techniques studied in Sections 121 and 01 

If we assume given an CTin^d) algorithm for multiplying two n by n polynomial matrices 
of degree d, combining the reductions of Section 0] with the cost estimates of Sections |2 
and 01 then yields 0~(n w <i) solutions to all our problems under consideration. Of course, we 
could have introduced a cost function MM(n, d) for polynomial matrix multiplication and 
derived more precise complexity estimates for each of the problems, in terms of (functions 
of) MM(n,d). However, we prefer for this paper to stick to the more readable 0~(n w <i) 
bound, which already gives a good sense of the link with polynomial matrix multiplication. 

A first task remaining would be to relax the regularity assumptions made for inversion 
(the input should be generic and of dimensions a power of two, see Section 14. 1|) and for 
row-reduction (the input should be non-singular, see Section |4~^|) . But even these "generic" 
situations are enough for our purpose here of showing how to rely on minimal bases and 
matrix fraction expansions/reconstructions. 

Also, recently, other problems on polynomial matrices than those treated in this paper 
have been shown to have about the same complexity as polynomial matrix multiplication. 
An example is the problem of computing the Smith normal form and thus also the de- 
terminant, whose solution in ^H] gives us Theorem 13.11 However — and this is the second 
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task remaining — , the list of problems that can be solved in about the same number of 
operations as for polynomial matrix multiplication still has to be augmented. The question 
is particularly interesting for the problem of computing the characteristic polynomial and 
the Frobenius normal form, for which the best known solutions fUEj have cost (J~(n 2 ' 7 d) 
still greater than 0~{n UJ d). 

Notation and basic reminders. Here and hereafter log denotes the logarithm in base two 
and I n the n by n identity matrix. For a matrix A over K[x], we denote its value at x = 
by A(0). For d G N and a matrix F over K[[x]], F = mod x d means that each entry of 
F is a multiple of x d , and F mod x d means that we truncate F into a polynomial matrix 
where only powers in x strictly less than d appear. By size of a polynomial matrix over 
K[x] we mean the number of elements of K that are necessary to represent it. For example, 
M G K[x] nxm of degree d has size at most nm(d + 1) = 0(nmd). A polynomial matrix 
is said to be non-singular when it is square and when its determinant is a non identically 
zero polynomial. Two matrices A,R G K[x] nxra are unimodularly left equivalent when there 
exists U G K[a;] nxn such that det U is a non-zero constant — that is, U is unimodular — and 
when UA = R. 

2 Minimal approximant bases and minimal nullspace bases 

Our solutions for solving a class of polynomial matrix problems in about the same number 
of operations in K as for multiplying two polynomial matrices will fundamentally rely on 
computing minimal bases of K[x]-modules. The target complexity estimate 0~(n w d) is 
reached since the bases we use are small, with size 0(n 2 d) is most cases, and may be 
computed fast (see Theorem 12.21 below) . 

Definition 2.1 Let Ai be a K[x]-submodule of'K[x] n of dimension T>. A basis Ni, . . . , N?> G 
K[x] n of Ai with degrees 5\ < ■ ■ ■ < 5t> is called a minimal basis if any other basis of Ai 
with degrees d\ < ■ ■ • < dx> satisfies di > 8j for 1 < i < T>. The degrees 5i are called the 
minimal indices of Ai. 

In applications to multivariable systems, this definition follows the study of minimal 
polynomial bases of vector spaces in 0. The two important examples of such bases that 
we use in this paper are minimal approximant bases and minimal nullspace bases. The 
approximant bases are defined from a power series matrix F over K[[x]], the nullspace bases 
are computed as special approximant bases from a polynomial matrix F = A over K[x]. 

2.1 Minimal approximant bases 

Given a formal power series F G K[[x]] nxra and an order d G N, we take for Ai the set of 
all approximants for F of order d: 

Ai = {v G K[x] lxn : vF = mod x d }. 
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The minimal bases of M. are called minimal approximant bases for F of order d. Since 
A4 has dimension n, such bases form non-singular n x n polynomial matrices. These 
polynomial matrices further have degree up to d and their size is thus of the order of n 2 d. 

Theorem 2.2 [8 . Let F £ K[[a;]] nxn and d £ N. A minimal approximant basis for F of 
order d can be computed in (Tfa^d) operations in K. 

Our notion of minimal approximant bases is directly inspired by pQ with some adap- 
tations for fully reflecting the polynomial matrix point of view. The cost estimate of 
Theorem 12.21 is a matrix polynomial generalization of the recursive Knuth/Schonhage half- 
gcd algorithm for scalar polynomials [1311111 (see also §11.1]), that takes into account 
fast polynomial matrix multiplication. 

For a matrix A over K[x], we denote by di its ith row degree, that is, the highest degree 
of all the entries of the ith row of A. The row leading matrix of A is the constant matrix 
whose zth row consists of the coefficients of x di in the ith row of A. We recall from [TU1 
§6.3.2] that a full row rank A is row-reduced when its row leading matrix also has full 
rank. As a consequence of their minimality, minimal approximant bases have the following 
properties, which will be used in Section 12.21 when specializing approximants for power 
series matrices to approximants for polynomial matrices. 

Property 2.3 Let N be a minimal approximant basis for F of order d. Then, 
I. N is row-reduced; 

II. If v £ M. has degree at most d, then there is a unique u £ K[x] lxn such that v = uN. 
Furthermore, N has at least one row of degree at most d. 

Property I above is a consequence of the minimality of the basis [TTH Theorem 6.5-10]. 
Property II is the fact that the rows of N form a basis, together with the predictable degree 
property fUl Theorem 6.3-13]. 

2.2 Minimal nullspace bases 

Given a polynomial matrix A £ K[x] nxn of rank r, we now take 

M = {v £ K[x] lxn : vA = 0}. 

This is a K[x]-submodule of K[x] n of dimension n — r. Its bases are called minimal nullspace 
bases for A and form full rank (n — r) x n polynomial matrices. The minimal indices 
Si < • • • < S n - r (see Definition 12.1)1 are called the (left) Kronecker indices of A jl()l §6.5.4]. 
For any given degree threshold 5, we further define 

K = max{l < i < n — r : 5i < 5}. 

A corresponding family of n linearly independent vectors of degrees 8\ , ■ ■ ■ , 5 K is a family 
of minimal nullspace vectors of degree at most 5. The theorem below says that if F = A is 
a polynomial matrix then any minimal approximant basis for A of sufficiently high order 
actually contains a family of minimal nullspace vectors for A. 
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Theorem 2.4 Let A G K[x] nxn be of degree d. Let N be a minimal approximant basis for 
A of order 5 + d + 1. Then exactly k rows of N have degree at most 5; these rows are in 
the (left) nullspace of A and their degrees are the Kronecker indices 5\, . . . , 8 K . 

Proof. A row v of iV of degree bounded by S satisfies vA = mod x 5+d+1 , and using 
degvA < 5 + d, vA = 0. Let k be the number of such v's in N, from the definition of k 
and since N is non-singular, k < n. We now verify that k > n. We consider k linearly 
independent vectors v , of degrees 5i in the nullspace of A. From Property 12 .HI we have 
V\ = U\M and deduce that one row of N has degree bounded by b\. Now, if N has i — 1 
rows of degrees bounded by 5i, . . . , d> t -_i, then the same reasoning with Vi as for Vi shows 
that has a row of degree bounded by <5j, linearly independent with respect to the first 
% — 1 chosen ones. It follows that k > k rows of iV have degrees bounded by 8±, . . . , 5 K , 
and are in the nullspace of A. Hence k = k, and we conclude using Definition 12.11 and the 
minimality of the (Vs. ■ 

For some applications, a shifted degree may be introduced (see [2] and the references 
therein), and some aspects of Theorem 12.41 mav be generalized accordingly (see [3 Theo- 
rem 4.2] or Lemma 6.3]). 

Notice that if the Kronecker indices of A are all bounded by d then an entire minimal 
nullspace basis for A can already be computed fast: by Theorem 12 A\ it suffices to compute 
a minimal approximant basis for A of order 2d+l and, by Theorem 12.21 this computation 
can be done in time (Tfa^d). 

However, in the general case of unbalanced degrees, computing a nullspace basis fast is 
much less immediate and the method we shall give in Section l4~4l relies on the complexity 
result given below in Theorem 12 .51 The cost given here is the one of a randomized algorithm 
of the Las Vegas kind — always correct, probably fast. The algorithm outputs correct 
minimal vectors in time 0~(n w d) with good probability, say greater than 1/2, otherwise 
returns failure (a correct result will be obtained after repetition). 

Theorem 2.5 17]. Let A £ «[x]^ n+m ^ xn with m < n be of full column rank and degree 
bounded by d. If 5 G N satisfies 

5m = 0(nd), (1) 

then a family of minimal nullspace vectors of degree at most 5 can be computed by a ran- 
domized Las Vegas (certified) algorithm in 0~(n u d) operations in K. 

Note that the cost estimate CF^n^d) relies on the compromise (P) between the minimal 
nullspace vector degree bound S and the row dimension of matrix A. For example, when 
m = 1 one can compute a nullspace vector of degree as large as 0(nd), whereas when 
m = n one may compute up to n nullspace vectors of degree 0(d). Random values are 
introduced essentially through a random compression matrix P G K[x] nxm that allows to 
compute minimal vectors more efficiently using the matrix AP G Kfa;]*™" 1 "" 1 -^" 1 rather than 
directly from A G K[x] (n+m)xm (see [13 Proposition 5.4]). 
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3 Matrix fraction expansion and reconstruction 

Matrix fraction expansion and reconstruction will be key tools especially for the row re- 
duction and the nullspace problems. Fraction reconstruction is a useful tool in computer 
algebra (e.g. see §5.7] for scalar polynomials), that is directly connected to coprime 
factorization (see below, and fOJ Chapter 6] or ^3] and the references therein). 

For a polynomial matrix A that is non-singular at x = and a polynomial matrix B, 
the techniques of [TT)| Proposition 17] reduce the computation of parts of the power series 
expansion 

oo 

A~ l B = J2 F i xi 

to polynomial matrix multiplication. By parts of the expansion, we mean a given number 
of consecutive matrix coefficients F{. This is summarized in the following theorem. 

Theorem 3.1 Let A G K[x] nx " with A(0) non-singular, and B G K[x] nxm . Assume 

that A and B have degree bounded by d and let h G N be such that h = 0{nd). If 6 G N 
satisfies 

5m = 0{nd), (2) 

then the 5 coefficients F^, F^ + i, . . . , Fh + s-\ G K nxm of the expansion of A~ l B at x = can 
be computed in 0~(n' J (i) operations in K. 

Similarly to Theorem 12 .5[ the cost estimate 0~(n^d) relies on the compromise (J2J) be- 
tween approximation order 5 and the column dimension of matrix B. For instance, for a 
vector B = b G K[x] nxl and h = 0, one can expand up to order 0(nd), whereas with 
B = I n and h = 0, one gets the expansion of A -1 up to order 0(d). In Section [Q| we 
shall use this result with B = I n and h = (n — l)d + 1 in order to get a high-order slice of 
length 0(nd) of the expansion of A^ 1 . 

Notice also that the regularity assumption det^4(0) ^ in Theorem 13 .11 is not restric- 
tive. Indeed, it can be satisfied with high probability using random shifts, thus yielding 
randomized algorithms for any A(0). Typically, with a randomly chosen xq G K, we shift 
x in the input like x <— x + x to get a regular input at zero and, at the end of the com- 
putation, we shift x back like x <— x — xq to recover the result (see [13 El HTj). 

A rational matrix H G K(x) nxm is strictly proper if lrnx^oo H(x) = G K nxm . In most 
applications, difficulties arise when A~ x G K(x) nxn is not strictly proper. However, one 
can define another fraction that is always strictly proper and shares some invariants with 
A^ 1 . Before seeing this, we first need to recall some facts about greatest common divisors 
of two polynomial matrices. 

Definition 3.2 A (left) matrix gcd of A £ K[x] nxn and B G K[x] nxm is any full column 
rank polynomial matrix G such that [G 0)U — [A B] with U unimodular over K[x]. 
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Definition 13.21 is for instance from [TU1 Lemma 6.3-3]. If [A B] has full row rank then 
all the gcd's of A and B are non-singular and equivalent with respect to multiplication on 
the right by any unimodular matrix in K[x] nxn (see 1 Q] Lemma 6.3-4]). A non-singular 
A G K[x] nxn is said to be (left) coprime with B G K[x] nxm if any gcd of A and B is 
unimodular (the gcd may be chosen as being the identity matrix /„). Similar definitions 
hold for rights gcd's and right coprimeness. 

Theorem 3.3 [Hj. Let A G K[x] nxri of degree bounded by d, with A(0) non-singular. For 
A' 1 = Y.T=o F i xi and h> (n- l)d, let H G K(x) nxn be given by H = ^Zo^+i^- Then 
H = A -1 (AH) = (HA) A" 1 is strictly proper, and AH and HA are polynomial matrices 
that are respectively left and right coprime with A. 

Proof. Let B = AH. By definition of H we have I n = A(A~ 1 mod x ) + x h B which 
in is (17) on the left with B and T respectively set to I n and A. It follows that 
B is a polynomial matrix. On the other hand, H = A _1 B is strictly proper because 
A~ l B = x~ h A~ 1 - x^ h (A^ mod x h ) and h > (n — l)d > deg A* where A* is the adjoint 
matrix of A. For establishing coprimeness we use 



[A x h B) 



I„ (A 1 mod x h 
L 



1 n 



In -A 



[In 0], (3) 



and the fact that if G is a left gcd of A and B it satisfies 

[G 0]U = [A B] (4) 

with U unimodular. Identities and (@J) give that their exists a polynomial matrix V 
such that [G 0}V = [I n 0], hence a polynomial matrix W such that GW = I n . Since G 
is a polynomial matrix this implies that G is unimodular, and A and B are left coprime. 
With B = HA, one could show similary right coprimeness. ■ 

For our application in Section l4~3^ we will need only the first, say 6, coefficients of the 
expansion of H as in Theorem 13.31 These coefficients thus correspond to a slice of order h 
and length 5 of the expansion of A~ l and, to recover them, we shall use Theorem 13 . 1 1 wit h 
B = I n . 

Matrix power series expansion will be used in conjunction with matrix (irreducible) 
fraction reconstruction or, equivalently, ( coprime ) factorization. We show below that min- 
imal approximant bases are appropriate tools for solving these problems. 

Definition 3.4 A (left) factorization of degree 5 of a rational matrix H G K(x) nxn is a 
representation H = V~ l U with U and V two polynomial matrices of degree bounded by 5. 
This factorization is said to be coprime when U and V are (left) coprime. 

A similar definition holds on the right. Hence, given H G K(x) nxri , the reconstruction 
or factorization problem is to recover U and V over K[x] such that V~ l U = H. If H is 
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defined at x = and given by its formal expansion F G K[[x]]™ xri , this problem reduces to 
computing a suitable [U V] G K[x]" x2n such that 



[U V] 



'In 



0. 



Theorem 3.5 Let H G K(x) nxn be strictly proper, with expansion F G K[[x]] 7ixra at x = 0. 
Assume that H admits a right factorization of degree 5r and a left factorization of degree Si- 
Let N G K[x] 2nx2n be a minimal approximant basis for [— I n F T ] T of order Sl+Sr+1. Then 
exactly n rows of N have degree bounded by 5l; these rows form a matrix [U V] G K[x] nx2ri 
such that V~~ X XJ is a left coprime factorization of H , with V row-reduced. 

Proof. Let BA~ l be a right factorization of H of degree 5r and T -1 ^ be a left factoriza- 
tion of H degree 5l- Since [—I n F T ] T A = [—A T B T ] T , iV is also a minimal approximant 
basis of the latter matrix whose rank is n. Using [S T}[-A T B T f = 0, with the threshold 
5 = 5l we have k = n. (See before Theorem 12.41 for a definition of k.) Hence, applying 
Theorem 12.41 to [— A T B T ] T (augmented on the right with n zero columns) with 5 = 5l 
and d = 5r, we know that exactly n rows of N have degree bounded by 5l and are in the 
nullspace of [— A T B T ] T . We denote the corresponding matrix by [U V]. The matrix V 
is non-singular, for otherwise there would be a non-zero vector v such that vV = 0. This 
would imply vVB = vUA = 0, hence either vU = or wA = for w = vU ^ 0, and would 
contradict either that rank[C/ V] = n or that A is non-singular. Therefore, V~ X U is a left 
factorization of H. 

This factorization must further be left coprime. Indeed, non-coprimeness would imply 
that U and V have a non-trivial left gcd, that is, there exists a polynomial matrix G such 
that U = GU', V = GU' and deg(det G) > 0. Then [GU 1 GV ] would be a submatrix of 
the minimal approximant basis, which would contradict its irreducibility in Theorem 
6.5-10] by considering a zero of det G. In addition, the fact that [U V] as a submatrix of iV 
is row-reduced (see Property 12. 3|) . implies that V is row-reduced. Indeed, since H = V~ 1 U 
is strictly proper, the row degrees of U are strictly smaller than those of V ^Ql Lemma 
6.3-10], and the row leading matrix of [U V] has the form [0 L\ where L is the row leading 
matrix of V, which is then non-singular. ■ 

As an immediate consequence of Theorem 13.51 and Theorem 12 .2\ coprime factorizations 
can be computed fast when the input matrix fractions admit left and right factorizations 
of degree 0(d). This corollary, given below, will be applied in Section [Ql to the particular 
matrix fraction of Theorem 13.31 

Corollary 3.6 Let H G K(x) nxn be as in Theorem [Q with S L = O(d) and S R = O(d). 
Given the first Sl + Sr + 1 coefficients of the expansion of H at x = 0, one can compute a 
left coprime factorization of H in 0~(n w d) operations in K. 
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4 Applications 

In this section, we show how the techniques presented in Sections El and El can be used to 
solve the following problems asymptotically fast: 

• \r\v n ^- given a non-singular A G K[x] nxri of degree d, compute A^ 1 . 

• Det n d : given A G K[x] nxn of degree d, compute det A. 

• RowRedn^: given A G K[x] nxn of degree d, compute a row-reduced form of A. 

• Nullspace„ d : given A G K[a;] nxn of degree d, compute the rank r of A and a full rank 
N G K[x]( n - r ) x " such that NA = 0. 

• Factor^: given a right factorization of degree d of H G K(x) nxn , compute a left 
factorization of H. 

Our approach here is to reduce each of the above five problems to (collections of) the 
problems below, for which (T(n w d) solutions are known: 

• MatMuU^: given A, B E K[x] nxn of degree d, compute the product AB. 

for solutions in time 0~{n u d) see |3]/, J^. 

• PartialNullSpace m s : given 5 = 0(nd/m) with n,d fixed, and given A £ K[x]^ n+m ^ xn 
of degree d, compute the minimal nullspace vectors of A of degree at most S. 

solved in time 0~(n u d) by Theorem \2.bl . 

• MatFracExp m s : given 5 = 0(nd/m) with n,d,h fixed such that h = 0(nd), and 
given A G K[x] nxn ,B e K[x] nxm of degree d with A(0) non-singular, compute the 5 
coefficients F h , F h+1 , . . . , F h+ g-i of the expansion of A~ X B at x — 0. 

> solved in time (Tiji^d) by Theorem \3.1\ 

• MatFracRec n(i : given 5l,5r = 0(d) and the first 5l + 5r + 1 coefficients of the 
expansion at x = of H G K(x) nxn as in Theorem 13 .5\ compute a left coprime 
factorization of H with row-reduced denominator. 

solved in time 0~(n"d) by Corollarv \3.b\ 

Assuming that n is a power of two and given a problem P n ^ or P mj $ such as any of those 
just introduced, we define the collections of problems we shall rely on as 

p n,d ■= ( solve °( 2i ) problems Pn/ 2%2 -d} < i<logn • (5) 

Such collections can be solved at about the same cost as polynomial matrix multiplication, 
as shown below. Here subscripts n, d and m, 5 should be added to P and P* depending on 
the underlying problem. 
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Lemma 4.1 For all P G {MatMul, PartialNullSpace, MatFracExp, MatFracRec}, one can 

solve P* in (T(n w d) operations in K. 

Proof. This an immediate consequence of (J^J) and of the bound O~(n a; o?) on the cost of 
each of these four problems. ■ 

4.1 Polynomial matrix inversion (lnv n d) 

Given A G K[x]" x ™ non-singular of degree d, the problem is to compute A^ 1 G K(x) nxn . 

Assuming that A is generic and that n is a power of two, we recall from ||| how 
lnv n rf reduces to PartialNullSpace* d plus some polynomial matrix multiplications. The 
algorithm in p.75] essentially consists in computing in logn steps a non-singular matrix 
U G K[x] nxn and a diagonal matrix B G K[x] nxn such that 

UA = B. (6) 

The inverse of A is then recovered as A^ 1 = B~ X \J ' . The first step is as follows. Let 
A = [A L A R ] where A L ,A R G K[x] nxn/2 and let N,Ne K[x] n / 2xn be minimal nullspace 
bases for, respectively, Al,Ar. This gives the first block-elimination step towards the 
diagonalization of A: 

A = [ A L A R ] -> NA = 

When A is generic of degree d, it turns out that all the minimal indices of both N_ and N 
are equal to d [HI Fact 1] and that NA^ and N An are n/2 x n/2 polynomial matrices of 
degree exactly 2d on which we iterate. 

We show in H that the property "dimension x degree = nd" generically carries from 
one iteration to the other: at step i, starting from 2 l ~ l blocks of dimensions (n/2 l ~ r ) x 

(n/2 t ~ 1 ) and degree 2 t ~ 1 d, we compute 2*" 1 pairs (N^\ N { 3 ) of minimal nullspace bases 
of dimensions (n/2*) x (n/2 l ~ 1 ) and whose minimal indices are all equal to 2*~ 1 <i. Let 
(U, B) = (/„, A) before the first step. Step i also requires to update the matrix transform 
as U <— diag[iV^]j x U and the right hand side as B <— di&g(N-^)j x B. Because of the 
special block-structure of the polynomial matrices involved, it can be shown that these 
updates reduce to solving 0{2 2t ) problems MatMul^i-i^- 1 ^ 

Overall, the logn block-diagonalization steps thus reduce to PartialNullSpace* d and to 

{solve 0(2 21 ) problems MatMul n/2 i j2 i d } <. <logn . (8) 

By Lemma f4. II and (jHJ), we therefore obtain a solution to lnv n ^ in (T(n 3 d) operations in K. 

Since by Cramer's rule each entry of A" 1 has the form p/(det A) where p G K[x] may 
have degree at large as (n — l)d, the size of A -1 is of the order of n 3 d. The above inversion 
algorithm, defined for A generic and n a power of two, is therefore nearly optimal. 



N 
N 



[A L A R ] 



NA, 



NA 



(7) 
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4.2 Determinant computation (Det n d) 



Given A £ K[s]" xri 0/ degree d, the problem is to compute det A £ K[x]. 

We assume here that A is generic with n is a power of two, and we use the inversion 
algorithm of Section |4~T1 It has been shown in 8 a that the diagonal entries of the diagonal 
matrix B in (jSJ) are constant multiples of det A. Since det A(0) is generically non-zero, we 
have 

, , det ,4(0), 

det A = — — -V^; i for all 1 < % < n. 

bi,i{0) M ~ ~ 

The problem Det n d thus reduces essentially to computing the determinant of the constant 
matrix A(0) and to the computation of, say, b\^. It is well-known that over K computing 
the determinant reduces to matrix multiplication jH Section 16.4] (that is, Det nj o reduces 
to MatMul ni o using our notations). Concerning bi t x, we perform logn steps as for inversion 
but, since &11 is the upper-left corner of B, we use instead of (0) the simpler step 

A= [A L A R ] -> NA L . (9) 

As in (JJJ), N is a minimal nullspace basis for Ar. Step i now consists in computing a 
single minimal nullspace basis of dimensions (n/2*) x (n/2 J_1 ) and minimal indices 2 l ~ 1 d, 
and then in multiplying this basis with the left half of an n/2 t by n/2 l block of degree 
2 l ~ 1 d, as in (jUJ). Hence, computing by performing these logn steps reduces to solving 
PartialNullSpace* d and MatMul* d . By Lemma l4~Tl this gives a solution to Det n)(i in 0~(n ul d) 
operations in K. 

Notice that when A is not generic or when n is not a power of two, a Las Vegas (Tfo^d) 
solution to Det n d can be obtained using the Smith normal form algorithm in [TK] . 



4.3 Row reduction (RowRed n d) 

Given A e K[a;] nxn of degree d, the problem is to compute R £ K[x] nxn that is row-reduced 
and unimodularly left equivalent to A. 

We assume here that A(0) is non-singular. Recall from Section 12.11 and ^01 §6.3.2] 
that R = A is a row-reduced form of A when R is row-reduced and R = UA for some 
unimodular polynomial matrix U. The solution in [H] works by expansion/reconstruction 
of the matrix fraction H as in Theorem 13.31 with h = (n — l)d+l. 

First, we expand H up to order 2d+ 1. This is done by solving MatFracExp n 2d+1 once, 
taking B = I n and h = (n — l)d + 1 = 0(nd). From Theorem 13.31 we know that H is 
a strictly proper matrix fraction which admits left and right factorizations A^ 1 {AH) and 
(HA)A~ 1 . Strict properness further implies that the degrees of both AH and HA must 
be less than the degree of A [TUl Lemma 6.3-10], and are thus bounded by d as well. 
Therefore, these left and right factorizations of H are factorizations of degree d and, using 
Theorem 13.51 we can reconstruct H from its expansion up to order 2d + 1 as H = R~ 1 S. 
This reconstruction corresponds to solving problem MatFracRec n) d once. On one hand, we 
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know by Theorem 13.51 that R is row-reduced. On the other hand, A~ l (AH) and R~ 1 S are 
coprime factorizations of the same fraction, which implies that there exists a unimodular U 
such that UA = R |10| Theorem 6.5-4]. It follows that R is indeed a row- reduced form of 
A. By Lemma [4. 1\ this reduction to MatFracExp n 2d+1 and MatFracRec n) d gives a solution 
to RowRed n; d in 0~{n ul d) operations in K. 

4.4 Small nullspace computation (Nullspace n d ) 

Given A G K[a;] nxn of degree d, the problem is to compute the rank r of A, and N G 

Kj^j (n-r)xn Q j ran fc n — r suc h that NA = 0. 

As already seen, a solution in the restrictive (e.g. generic) case when all minimal vectors 
have degrees in 0(d) is provided by a solution to PartialNullSpace n d . In the general case 
the row degrees in a nullspace basis of A may be unbalanced, they range between and 
nd fTf\ Theorem 3.3]. Previously known methods, whose cost is essentially driven by the 
highest Kronecker index, do not seem to allow the target complexity estimate 0~(n u3 d) (see 
for instance [13 Section 2]). 

Our solution in [Tj\ first reduces the general nullspace problem to the full column rank 
case via randomization. This consists in evaluating the rank r of A at a random x = xq, 
then in compressing A to a full column rank matrix. We also derive a particular strategy 
when n ^> r. Consequently, for a simplified explanation here, we now assume that A has 
full column rank n and dimensions (n + m) x n with m = 0(n). 

The algorithm then works in % steps with 1 < i < logn. At step % we compute a set 
of about mjT nullspace vectors of degrees less that S = 2 l d. These vectors are obtained 
from logn solutions to PartialNullSpace, m s for nullspace vectors of bounded degree 5 = 2 l d, 
and involving matrices of decreasing dimensions n + mjT . Hence we essentially have a 
reduction to PartialNullSpace^ s . We may point out that the proof of Theorem 12.51 for the 
cost of the partial nullspace itself relies on solutions to MatFracExp m s , and MatFracRec m ^. 
Nullspace vectors are computed using a matrix fraction expansion /reconstruction scheme. 

The appropriate instances for PartialNullSpace m/ / 2 \2M5 1 — * — logn, are built as subma- 
trices of the input matrix A. Our choices for these submatrices ensure the linear indepen- 
dency of the successive computed sets of nullspace vectors. The algorithm hence outputs a 
union of a logarithmic number of sets of linearly independent nullspace vectors. Each set, 
corresponding to an instance of PartialNullSpace m / 2 , is a family of minimal vectors for 
a submatrix of A. The minimality is not preserved in general with respect to A, however 
we prove that small degree vectors are obtained fT7\ Proposition 7.1]. 

This reduction of NullSpace n d to PartialNullSpace^ s and to MatMul* d for additional 
matrix multiplications establishes that a solution matrix iV such that NA = can be 
computed in 0~(n w d) operations in K by a randomized Las Vegas (certified) algorithm. 

4.5 Factorization (Factor^) 

Given a right factorization BA~ l of degree d of H G K(x) nxn , the problem is to compute 
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polynomial matrices U and V such that V 1 U = H . 



Corollary 13 .6[ together with the expansion of H = BA" 1 , provides a solution to 
FracMatRec n)( i if H admits factorizations of degree d on both sides. The solution of the 
general case, we mean for an arbitrary left side factorization, induces several difficulties for 
dealing with unbalanced row degrees. These difficulties are bypassed using the techniques 
of Section IP1 

By considering the polynomial matrix [— A T B T ] and solving Nullspace 2ni d we get U 
and V such that 

" —A 
B 



[U V] 



0. 



Arguments similar to those used in the proof of Theorem 13.51 lead to the fact that V is 
non-singular. Hence a solution V~ X U to the factorization problem is computed in 0~(n w d) 
operations in K. Note that since a solution to Nullspace2 n ,d rnay not be minimal, the 
factorization V^U may not be coprime. 
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